Models of universe with a polytropic equation of state: 
II. The late universe 
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We construct models of universe with a generalized equation of state p = (ap + kp 1+1 ^ n )c 2 
having a linear component and a polytropic component. The linear equation of state p = ape 2 
describes radiation (a = 1/3), pressureless matter (a = 0), stiff matter (a = 1), and vacuum 
energy (a = —1). The polytropic equation of state p = kp 1+1 ' n c 2 may be due to Bose-Einstein 
condensates with repulsive (k > 0) or attractive (k < 0) self-interaction, or have another origin. In 
this paper, we consider negative indices n < 0. In that case, the polytropic component dominates 
in the late universe where the density is low. For a = 0, n = —1 and k — —pA, we obtain a model of 
late universe describing the transition from the matter era to the dark energy era. Coincidentally, 
we live close to the transition between these two periods, corresponding to a 2 = 8.95 10 25 m and 
ta = 2.97 10 17 s. The universe exists eternally in the future and undergoes an inflationary expansion 
with the cosmological density pa = 7.02 1CP 24 g/m 3 on a timescale tA = 1.46 10 18 s. For a = 0, 
7i = —1 and k — pA, we obtain a model of cyclic universe appearing and disappearing periodically. 
If we were living in this universe, it would disappear in about 2.38 billion years. We make the 
connection between the early and the late universe and propose a simple equation describing the 
whole evolution of the universe. This leads to a model of universe that is eternal in past and future 
without singularity (aioniotic universe). It generalizes the ACDM model by removing the primordial 
singularity (Big Bang). This model exhibits a nice "symmetry" between an early and late phase 
of inflation, the cosmological constant in the late universe playing the same role as the Planck 
constant in the early universe. We interpret the cosmological constant as a fundamental constant 
of nature describing the "cosmophysics" just like the Planck constant describes the microphysics. 
The Planck density and the cosmological density represent fundamental upper and lower bounds 
differing by 122 orders of magnitude. The cosmological constant "problem" may be a false problem. 
We determine the potential of the scalar field (quintessence, tachyon field) corresponding to the 
generalized equation of state p = (ap + kp 1+1/ln )c 2 . We also propose a unification of pre-radiation, 
radiation and dark energy through the quadratic equation of state p/c 2 — —4p 2 /3pp + p/3 — 4/9a/3. 
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I. INTRODUCTION 



According to contemporary cosmology, the present en- 
ergy content of the universe is composed of approxi- 
mately 5% ordinary matter, 20% dark matter and 75% 
dark energy The expansion of the universe began 
in a tremendous inflationary burst driven by the vacuum 
energy with the Planck density pp — 5.16 10" g/m 3 . Be- 
tween 10 -35 and 10 -33 seconds after the beginning (Big 
Bang), the universe expanded by a factor 10 30 @, In- 
flation does not offer any explanation for the time before 
that "beginning" . The universe then entered in the radi- 
ation era and, when the temperature cooled down below 
approximately 10 3 K, in the matter era At present, it 
undergoes an accelerated expansion Q presumably due 
to the cosmological constant or to some form of dark en- 
ergy with negative pressure violating the strong energy 
condition Q . This corresponds to a second period of in- 
flation, which is different from the first since it is driven 
by the cosmological energy density pa — 7.02 10~ 24 g/m 3 
instead of the Planck density. Despite the success of this 
model, the nature of dark matter, dark energy, and of 
the very early universe (pre-radiation era) remains very 
mysterious and leads to many speculations. 

The phase of inflation in the early universe is usu- 
ally described by some hypothetical scalar field <f> with 
its origin in quantum fluctuations of vacuum Q. This 



leads to an equation of state p = —pc 2 , implying a con- 
stant energy density, called the vacuum energy. This 
energy density is usually identified with the Planck den- 
sity pp. As a result of the vacuum energy, the universe 
expands exponentially rapidly on a timescale of the or- 
der of the Planck time tp = 5.39 10 _44 s (early inflation). 
This phase of inflation is followed by the radiation era de- 
scribed by an equation of state p = pc 2 /3. In our previ- 
ous paper (Paper I) , we showed that a unified description 
of the pre-radiation (p = pp) and radiation (p oc a~ 4 ) 
eras, connected by a phase of exponential inflation, could 
be obtained from a single equation of state of the form 
p = (l/3)/o(l — 4p/pp)c 2 . This provides a non-singular 
model of the early universe. 

The phase of acceleration in the present universe is 
usually ascribed to the cosmological constant A which is 
equivalent to a constant energy density p\ = A/(8nG) 
called dark energy. This can be modeled by an equation 
of state p = —pc 2 , implying a constant energy density 
identified with the cosmological density p\. As a result 
of the dark energy, the universe expands exponentially 
rapidly on a timescale of the order of the cosmological 
time t\ = 1.4610 18 s (de Sitter solution). This leads 
to a second phase of inflation (late inflation). Inspired 
by the analogy with the early inflation, some authors 
have represented the dark energy by a scalar field called 
quintessence 0. As an alternative to the quintessence, 
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other authors [8j] have proposed to model the accelera- 
tion of the universe by an exotic fluid with an equation 
of state of the form p = —A/p called the Chaplygin gas 
(see for a complete list of references). At late times, 
this equation of state leads to a constant energy density 
implying an exponential growth of the scale factor that 
is similar to the effect of the cosmological constant. At 
earlier times, this equation of state returns the results 
of the cold dark matter model. Therefore, it provides 
a unification of dark matter (p cx a~ 3 ) and dark energy 
(p = pa) in the late universe. Furthermore, it gives a 
real velocity of sound which is non-trivial for fluids with 
negative pressure. Some generalizations of this equation 
of state have been considered in the form p = —A/p a 
with a > — 1 0, H[ . As mentioned in Q , the Chaplygin 
gas has some connection with string theory and can be 
obtained from the Nambu-Goto action for d-branes mov- 
ing in a (d + 2)-dimcnsional spacetime in the light-cone 
paramctrization. Furthermore, it is the only fluid which, 
up to now, admits a supersymmetric generalization. 

From a theoretical point of view, it is desirable to study 
models of universe with a generalized equation of state 
p = (ap + fcp 1+1 /")c 2 having a standard linear compo- 
nent and a polytropic component. The linear equation 
of state p = ape 2 describes radiation (a = 1/3), pres- 
sureless matter (a = 0), stiff matter (a = 1), and vac- 
uum energy (a = —1). The polytropic equation of state 
p = kp 1 c 2 with 7 = 1 + 1 /n may be due to Bose-Einstein 
condensates with repulsive (k > 0) or attractive (k < 0) 
self- interaction Q, or have another origin. When n > 0, 
the polytropic component dominates the linear compo- 
nent when the density is high. These models, studied in 
Paper I, describe the early universe. Conversely, when 
n < 0, the polytropic component dominates the linear 
component when the density is low. These models, stud- 
ied in the present paper, describe the late universe. As 
we shall see, these two studies are strikingly symmetric. 
Interestingly, this symmetry seems to reflect the true evo- 
lution of the universe. 

In this paper, we propose an exhaustive study of the 
equation of state p = (ap+kp 1+1 / n )c 2 with n < 0. When 
k > 0, the universe exhibits a future peculiarity: Its den- 
sity vanishes in infinite time (n < —2), or periodically in 
time (n > —2), while its radius tends to a constant value 
(the universe is singular when n > — 1 since the pressure 
becomes infinite in finite time). For a = 0, n = — 1 and 
k = pa, we obtain a model of cyclic universe, equivalent 
to the anti-ACDM model, in which the universe "disap- 
pears" (it becomes empty) and "reappears" periodically. 
According to this model, the universe would disappear 
in about 2.38 billion years. When k < 0, the universe 
exists for all times in the future and there is no singu- 
larity. For a = 0, n = —1 and k = —pA, we obtain a 
model of late universe, equivalent to the ACDM model, 
describing in a unified manner the transition from the 
matter era to the dark energy era. Coincidcntally, we 
live close to the transition between these two periods, 
corresponding to a,2 = 8.95 10 25 m and ti = 2.9710 17 s. 



This universe exists eternally in the future and undergoes 
an inflationary expansion with the cosmological density 
PA = 7.02 lCT 24 g/m 3 on a timescale i A = 1.46 10 18 s. 

The paper is organized as follows. In Sec. HH we recall 
the basic equations of cosmology. In Sees. IIIII and IIV1 
we study the generalized equation of state p = (ap + 



kp 



l + l/n\ 



for any value of the parameters a, k and n < 



0. In Sec. [V] we consider the case a = 0, n = — 1, and 
k = —pA, providing a model of non-singular inflationary 
universe describing the transition from the matter era 
to the dark energy era (ACDM model). In Sec. IVI1 we 
consider the case a = 0, n = — 1, and k = p\, leading to 
a model of peculiar cyclic universe (anti-ACDM model) . 
In Sec. IVIIl we discuss the connection between the early 
and the late universe, and propose a simple equation [sec 
Eq. (|97p] describing the whole evolution of the universe. 
This leads to a model of universe that is eternal in past 
and future without singularity (aioniotic universe). This 
model exhibits a nice "symmetry" between an early and 
late phase of inflation, the cosmological constant in the 
late universe playing the same role as the Planck constant 
in the early universe. This model refines the standard 
ACDM model by removing the primordial singularity. In 
Sec. IVIII1 we determine the potential of the effective 
scalar field (quintessence, tachyon field) corresponding to 
the generalized equation of state p = [ap + kp 1+1 ^ n )c 2 . 



II. BASIC EQUATIONS OF COSMOLOGY 

In a space with uniform curvature, the line element 
is given by the Friedmann-Lemaitre-Roberston- Walker 
(FLRW) metric 
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+ r 2 (d0 2 + sin 2 d(t> 2 ) 
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where a(t) represents the radius of curvature of the 3- 
dimensional space, or the scale factor. By an abuse of 
language, we shall call it the "radius of the universe". 
On the other hand, k determines the curvature of space. 
The universe can be closed (k > 0), flat (k = 0), or open 
(k < 0). 

If the universe is isotropic and homogeneous at all 
points in conformity with the line element ([1]), and con- 
tains a uniform perfect fluid of energy density e(t) = 
p(t)c 2 and isotropic pressure p(t), the Einstein equations 
can be written as 
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where we have introduced the Hubble parameter H = 
a/ a and accounted for a possibly non-zero cosmologi- 
cal constant A. These are the well-known Fricdmann 
equations describing a non-static universe Q. Among 
these three equations, only two are independent. The 
first equation can be viewed as an "equation of continu- 
ity". For a given barotropic equation of state p = p(p), 
it determines the relation between the density and the 
scale factor. Then, the temporal evolution of the scale 
factor is given by Eq. 

We will also need the thermodynamical equation 



§L - 1 

dT ~ T 



(5) 



which can be derived from the first principle of thermo- 
dynamics [H. For a given barotropic equation of state 
P — P(p)> this equation can be integrated to obtain the 
relation T = T(p) between the temperature and the den- 
sity. It can be shown Q that the Friedmann equations 
conserve the entropy of the universe 
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In this paper, we consider a flat universe (k = 0) in 
agreement with the observations of the cosmic microwave 
background (CMB) On the other hand, we set A = 
because the effect of the cosmological constant will be 
taken into account in the equation of state. The Fried- 
mann equations reduce to 



dp . „a 



AttG 



0. 



(7) 
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where k is the polytropic constant and 7 = l + l/nis 
the polytropic index. Concerning the linear equation of 
state, we assume — 1 < a < 1 (the case a = — 1 is treated 
specifically in Appendix [X| . This equation of state de- 
scribes radiation (a = 1/3), pressureless matter (a = 0) 
and vacuum energy (a = —1). Concerning the poly- 
tropic equation of state, we remain very general. Such 
an equation of state may correspond to self-gravitating 
BECs with repulsive (k > 0) or attractive (k < 0) self- 
interaction, but it may have another origin. 

For an equation of state of the form (fTTJ]) with posi- 
tive index n > 0, the polytropic component dominates 
the linear component when the density is high. These 
models, studied in paper I, describe the early universe. 
Conversely, when n < 0, the polytropic component dom- 
inates the linear component when the density is low. 
These models, studied in the present paper, describe the 
late universe. We shall see that these two situations are 
strikingly similar. They reveal a form of "symmetry" be- 
tween the early and the late universe. In this paper, as 
in Paper I, we assume that a + 1 + kp x l n > 0. This 
corresponds to the "normal" case where the density de- 
creases with the radius (w > —1). The opposite case 
q + 1 + kp 1 /' 1 < 0, leading to a "phantom universe" 
where the density increases with the radius (w < — 1), 
is treated in Paper III. 

The purpose of our series of papers is to make an ex- 
haustive study of all the possible cases (a,n,k), even if 
some of them seem to be unphysical or in conflict with 
the known properties of our universe. Actually, by re- 
jecting the cases leading to past or future singularities 
(or peculiarities), we will be led naturally to the "good" 
model (see Sec. NTTPl . 

Some of the equations established in Paper I remain 
valid when n < 0. To avoid repetitions, we shall directly 
refer to these formulae. For example, Eq. (1-60) refers to 
Eq. (60) of Paper I. 
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A. The density 



Introducing the equation of state parameter w = p/pc 2 , 
we see from Eq. |(5J) that the universe is decelerating if 
w > —1/3 (strong energy condition) and accelerating if 
w < —1/3. On the other hand, according to Eq. ([7]), the 
density decreases with the scale factor if w > —1 (null 
dominant energy condition) and increases with the scale 
factor if w < — 1 . 



III. GENERALIZED EQUATION OF STATE 
WITH NEGATIVE INDEX 

We consider a generalized equation of state of the form 

p=(ap + kp 1+1 / n )c 2 . (10) 

This is the sum of a standard linear equation of state 
p = ape 2 and a polytropic equation of state p = kp^c 2 , 



For the equation of state ([7]), assuming w > — 1, the 
Fricdmann equation ^ can be integrated into 



(11) 



[(a/a^t^/^l] 



where = [(«+ l)/|fc|]™ and a* is a constant of integra- 
tion. The upper sign corresponds to k > and the lower 
sign corresponds to k < 0. 

For k > 0, the density is defined only for a < a*. When 
a — >• 0, p/p* ~ (a*/a) 3 ( 1+a ) — > +oo corresponding to the 
linear equation of state. When a — > a* , 



p* 



3(1 + a) 
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(1 — a/a*) 
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For the behavior of the pressure, we must consider dif- 
ferent cases. When a — > 0, p — > +oo for a > 0, p — > — oo 
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FIG. 1. Density as a function of the scale factor for k > 
and k < 0. We have taken n = — 1/2 and a = 0. 



Some curves giving the evolution of the temperature T 
as a function of the scale factor a are plotted in Fig. [2] 
for fc > and fc < 0. 




for a < 0, p — > +oo for (a = 0,n < —1), p — >• for 
(a = 0, n > —1), and p = fee 2 for (a = 0,n = —1). When 
a — > a*, p — > +oo for n > — 1, p tends to a finite value 
for n = — 1, and p — > for n < — 1. 

For fc < 0, the density is defined for all a. When a — > 0, 
p ~ p*( a */ a ) 3 ^ 1+c ^ — > +oo corresponding to the linear 
equation of state. When a — > +oo, the density tends to 
a finite value p*. For the behavior of the pressure, we 
must consider different cases. When a — >• 0, p — s- +oo for 
a > 0, p — >• — oo for a < 0, p — > — oo for (a = 0, n < —1), 
p — > for (a = 0, n > —1), and p = fee 2 for (a = 0, n = 
— 1). When a — > +oo, p — > — p*c 2 . 

Some curves giving the evolution of the density p as 
a function of the scale factor a are plotted in Fig. [1] for 
k > and fc < 0. 



B. The temperature 

For the equation of state (fTU|) . the evolution of the 
temperature is given by Eqs. (F60) and (1-61). Let us 
consider some asymptotic limits. 

When a <C a*, Eq. (F61) reduces to Eq. (F62) re- 
turning the relation valid for a linear equation of state. 
When a 0, T -> +oo for a > 0, T -> for a < 0, and 
T^n for a = 0. 

When fc > and a —> a* . 



-| Tl + l 



3(1 
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(13) 



(1 - a/a*) n+1 

Tl for n = — 1 and 



Therefore, T -> for n < -1, T 
T -> +oo for n > — 1. 

When fc < and a — > +oo, T/T* is given by Eq. (I- 
64). Therefore, T -> for a + n + 1 > 0, T -> T* for 
a + n + 1 = and T -> +oo for a + n + 1 < 0. 

The extremum of temperature (when it exists) is lo- 
cated at the point defined by Eqs. (I-65)-(I-67). Finally, 
the entropy © is given by Eq. (1-68). 



FIG. 2. Temperature as a function of the scale factor for 
k > and k < 0. We have taken n = — 1/2 and a = 0. 



C. The equation of state parameter w(t) 



We can rewrite the equation of state (|T0|) as p = 
w(t)pc 2 where w(t) is given by Eq. (1-69). The point 
(a~w , Pw , T w ) where the pressure vanishes is defined by 
Eqs. (1-70) and (1-71). 

For fc > 0, w — > a when a — > 0, and w — > +oo when 
a — > a*. For a > 0, the pressure is always positive (w > 
0). For a < 0, the pressure is negative when a < a w and 
positive when a w < a < a*. 

For fc < 0, w — > a when a — > and u> — > — 1 when 
a — > +oo. For a < 0, the pressure is always negative 
(w < 0). For a > 0, the pressure is positive when a < a w 
and negative when a > a w . 
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FIG. 3. The parameters u;, g and c? a /c? as a function of the 
scale factor a for k > 0. We have taken n = — 1/2 and a — 0. 
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FIG. 4. The parameters w, q and c? a /c? as a function of the 
scale factor a for k < 0. We have taken n — — 1/2 and a = 0. 



Some curves giving the evolution of w as a function of 
the scale factor a are plotted in Figs. [3] and 0] for k > 
and fc < 0. 



D. The velocity of sound 

For the equation of state (fTU|) , the velocity of sound is 
given by Eq. (1-72). The velocity of sound vanishes at 
the point (a e ,p e7 T e ) defined by Eqs. (I-65)-(I-67) where 
the temperature is extremum. At that point the pressure 
is extremum with value p e given by Eq. (1-73). The case 
c 2 < corresponds to an imaginary velocity of sound. 
The point (a s , p s ,T s ) where the velocity of sound is equal 
to the speed of light is defined by Eqs. (1-74)- (1-76). 
Different cases have to be considered. 

We first assume k > 0. When a — > 0, (c s /c) 2 — > a; 
when a — > a*, (c s /c) 2 — > +oo for n < —1 and (c s /c) 2 — > 
—oo for n > — 1. For n < —1 and a > 0, c 2 s is always 
positive. For n < — 1 and a < 0, c 2 is negative when 
a < a e and positive when a e < a < a*. The velocity of 
sound is smaller than the speed of light when a < a s and 
larger when a s < a < a*. For n > — 1 and a < 0, c 2 in 
always negative. For n > — 1 and a > 0, c 2 is positive 
when a < a e and negative when a e < a < a*. For n > — 1 
and a < 1, the velocity of sound is always smaller than 
the speed of light. 

We now assume k < 0. When a —> 0, (c s /c) 2 —> a; 
when a — > +oo, (c s /c) 2 — > —(a+n+l)/n. For a+n+l > 
and a > 0, c 2 is always positive. For a + n + 1 > 
and a < 0, c 2 is negative when a < a e and positive when 
a > a e . For a + n + 1 < and a < 0, c 2 in always 
negative. For a + n + 1 < and a > 0, c 2 is positive 
when a < a e and negative when a > a e . For a < 1 and 
a + 2n + 1 < 0, the velocity of sound is always smaller 
than the speed of light. For a < 1 and a + 2n + 1 > 0, 
the velocity of sound is smaller than the speed of light 
when a < a s and larger when a > a s . 

Some curves giving the evolution of (c s /c) 2 as a func- 



tion of the scale factor a are plotted in Figs. [3] and |4] for 
k > and k < 0. 



IV. EVOLUTION OF THE SCALE FACTOR 
A. The deceleration parameter 

The deceleration parameter is defined by Eq. (1-77). 
In a flat universe, we obtain the relation of Eq. (1-78). 
The universe is decelerating when q > and accelerating 
when q < 0. For the equation of state ([TO]) , we get Eq. (I- 
79). The point (a c ,p c ,T c ) defined by Eqs. (1-80) and (I- 
81) corresponds to a possible inflexion point (q = a = 0) 
in the curve a(t). Different cases have to be considered. 

We first assume k > 0. When a -> 0, q -> (1 + 3a)/2; 
when a — > a*, g — > +oo. For a > —1/3, the universe is 
always decelerating (q > 0). For a < —1/3, the universe 
is accelerating when a < a c and decelerating when a c < 
a < a*. 

We now assume fc < 0. When a — > 0, g — > (1 + 3a)/2; 
when a — > +oo, — > — 1. For a < —1/3, the universe is 
always accelerating (q < 0). For a > —1/3, the universe 
is decelerating when a < a c and accelerating when a > 
a c . 

Some curves giving the evolution of q as a function of 
the scale factor are plotted in Figs. [3] and 2] for k > 
and k < 0. 



B. The differential equation 

The temporal evolution of the scale factor a(t) is de- 
termined by the Fricdmann equation Introducing 
the normalized radius R = a/a*, the density (|11[) can be 
written as Eq. (1-82) . Substituting this expression in Eq. 
dU, we obtain the differential equation 

a _ eKR 

[_R3(l+c0/n zp !]n/2 ' t X4 -> 

where if = (SttG/t^/S) 1 / 2 and e = ±1. In general, we 
shall select the sign e = +1 corresponding to an expand- 
ing universe (R > 0), except in the case of a cyclic uni- 
verse where both signs of e must be considered. The solu- 
tion of Eq. (|14[) can be written as in Eq. (1-84) or, after 
a change of variables x = R 3 ( 1+a ')/ n : as in Eq. (1-85). 

The solution of Eq. (fT4"]) can be expressed in terms 
of hypergeometric functions. Some simple analytical ex- 
pressions can be obtained for specific values of n. Ac- 
tually, we can have a good idea of the behavior of the 
solution of Eq. (|14[) by considering asymptotic limits 
(see below). The complete solution is represented in the 
figures by solving Eq. (|14[) numerically. 
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C. The case k > 

For k > 0, the universe starts at t = with a vanishing 
radius R — 0, an infinite density and an infinite pressure 
(Big Bang singularity). When t — > 0, we recover the 
solution corresponding to a linear equation of state 



R 



3(a + l) 



-,2/[3(l+a)] 



A'f 



[§(a + l)Ai]' 



(15) 



(16) 



The final fate of the universe depends on the value of n. 
Different cases must be considered. 

For n < —2, the radius tends to a finite value R = 1 in 
infinite time, while the density and the pressure decrease 
to zero. When t — > +oo, 



1-R 



-2 |_3(a + l) 



H/2 



2/(|n|-2) 



(17) 
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FIG. 5. Evolution of the radius R as a function of time for 
k > 0. We have represented the cases n < —2, n = —2, 
and n > —2 (specifically n = —3, n = —2, and n = —1/2). 
We have taken a — 0. The universe evolves towards a future 
peculiarity in which the density tends to zero while the radius 
tends to a finite value (see text for details). For n > — 1, the 
pressure is infinite at t = t*. For n < — 1, the velocity of 
sound becomes larger than the speed of light when R > R s — 
{{a + 2n + l)/[n(l - a)]} n/[3(1+Q)1 . 



/->* 



2 3(a + l)^ 



2M/(M-2) 



(18) 



The density decreases algebraically rapidly. 

For n = — 2, the radius tends to a finite value U = 1 in 
infinite time, while the density and the pressure decrease 
to zero. When t —> +oo. 



1 - R ~ Ae 



3(l+cQ 



Kt 



P* 



n 2 



■(l + a)A 



,-3(l+a)Kt 



(19) 



(20) 



The density decreases exponentially rapidly. 

For n > —2, the evolution of the universe is cyclic. 
Starting from R = at t — (Big Bang), the radius 
increases, reaches its maximum value R = 1 at t = t*, 
then decreases and reaches zero at t = 2t* (Big Crunch), 
increases again and so on (since a = oo at t = 2t*, it may 
not be allowed to extend the solution after 2t t ). The 
density is infinite when R = and vanishes when R = 1. 
When t-*U, 



1-R 



P* 



3(a + l) 



M/2 



2/(2-|n|) 



K{U - 1) 



\n\ 3(a< + l) 



2|n|/(2-|n|) 



(21) 



(22) 



The universe displays a future peculiarity at t = t* since 
its density vanishes (the universe "disappears" ) while its 
radius has a finite value R = 1. At 4 = t„, the pressure 
tends to zero for n < — 1, is constant for n = — 1, and 



is infinite for n > — 1. Therefore, for n > — 1 there is a 
future singularity (since o = oo at £ = £* , it may not be 
allowed to extend the solution after t* in that case). 

The evolution of the scale factor is represented in Fig. 
[5] in the different cases described above. Some simple 
analytical results can be obtained in particular cases. 

For n = — 1, using the identity 



1 



dx 



\/x - 1 x 



2 arctan [yx 



1 



we obtain 



- ^2/[3(l+a)] 



R = sin 

P_ 

P* 



-{l + a)Kt 



tan 2 [\{l + a)Kt] 



(23) 



(24) 



(25) 



This describes a peculiar cyclic universe (see Sec. I"VT|) . 
The radius reaches its maximum value R = 1 with a 
vanishing density p = for the first time at Xt* = 
7r/[3(l + a)]. We can explicitly check that Eq. JM]) has 
the asymptotic forms (fTS"]) and (|2T1) . 
For n = — 2, using the identity 



1 



dx l x 
In 



1 x 



1 



we obtain 



R = 



\(l+a)Kt 



2/[3(l+a)] 



P* 



(26) 



(27) 



(28) 
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This describes a universe reaching a peculiarity R = 1 
with p = in infinite time. Equation (|27[) has the asymp- 
totic forms (15) and (H) with A = 2/[3(l + a)]. 
For ?i = —1/2, we have the identity 

1 (to /r L /r, 

V2(a; 



(.t-1)V4 x 



= v^arctan 



1 



arctan 1 — \/2( 



1)1/4 
1)1/4 



In 



-1 + V2(a; - l) 1 / 4 - y/^T 
1 + V2(x - l)i/ 4 + a/e^T 



(29) 



which determines i(i?) using Eq. (1-85). For a = 0, this 
solution corresponds to the anti-Chaplygin gas p = A/ p. 



D. The case k < 

For fc < 0, the universe starts at t = with a vanishing 
radius R = 0, an infinite density, and an infinite pressure 
(Big Bang singularity). When t — > 0, we recover the 
solution (|15[) -(fT6 ) corresponding to a linear equation of 
state. When t —> +oo, the density tends to a constant 
value leading to an exponential growth of the radius 



R~ Ae 1 



(30) 



This late inflationary phase can mimic the effect of the 
cosmological constant. This was the original motivation 
of the Chaplygin gas model corresponding to a = and 
n = —1/2 The evolution of the scale factor is repre- 
sented in Fig. |6] Some simple analytical results can be 
obtained in particular cases. 




p* tanh 2 [§(1 + a)Kt] 



(33) 



We can explicitly check that Eq. (|52"|) has the asymptotic 
forms dHJ) and (3D) with A = i/2 2 /[ 3 (i+ Q )l. 
For n = — 2, using the identity 



1 dx 



x + 1 x 



In 



1 +x 



we obtain 



i? = 

/9* 



[(l+a)Jfi _ 1 



2/[3(1+q)] 



1 - e ~ 3 



(34) 



(35) 



(36) 



Eq. ([33]) has the asymptotic forms (|15p and ([3D) with 
A= 1. 

For n = —1/2, we have the identity 



(X + 1)V 4 X 



dx 

— = 2 arctan 



+ ln 



(1 + x) 1 / 4 
(l+x) 1 / 4 -! 



(37) 



which determines t(^) using Eq. (1-85). For a = 0, this 
returns the result of [8f for the Chaplygin gas p = —A/ p. 



V. A MODEL OF LATE INFLATIONARY 
UNIVERSE 

The generalized equation of state (flO) with n < and 
fc < can be used to describe the transition between a 
phase where the universe has a linear equation of state 
(p = ape 2 ) and a phase of late inflation where its density 
is constant (p = p*). Therefore, if we take a = 0, it 
describes the transition between the matter era and the 
dark energy era. This provides a model of late inflation- 
ary universe. For simplicity, we shall take n = — 1. In 
this way, we obtain a model of late universe that is "sym- 
metric" to the model of early universe studied in Paper 
I, corresponding to an index n = 1. 



FIG. 6. Evolution of the radius 7? as a function of time for 
k < 0. We have taken n = —1/2 and a = 0. This corresponds 
to the original Chaplygin gas @]. 



For n = —1, using the identity 



1 dx , / a/1 + a; - 1 
— = m 



y/x + 1 x 



we obtain 



R = sinh 



2/[3(l+a)] 



-(l+a)Jft 



(31) 



(32) 



A. The basic equations 

For a = 0, n = -1, and k < 0, Eqs. (JTD) and (|TT) 

become 



p = kc 



P = P* 



^] +1 

a 



(38) 



(39) 



where = |fc| and we have written 02 for a*. 
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When a > 02, the density is approximately constant: 
p ~ . This gives rise to a phase of late inflation. The 
value p„ defines a fundamental lower bound p m in for the 
density. We shall identify it with the dark energy density 
p\ = 7.02 1 - 24 g/m 3 (see Appendix [Bjl . Hence, we take 



P* — Pmin — PA- 



(40) 



This fixes k = —p\. This is how the dark energy (cos- 
mological constant) is taken into account in our model. 

When a<fl2, Eq. (f39|) reduces to p ~ PA^/a) 3 cor- 
responding to a pressureless universe (p = 0). This de- 
scribes the matter era. The conservation of p m o? implies 
p K al = p mfi al. Writing p mfi = fl m ,oPo and p A = Oa,oPo 
where il m ,o = ^b,o + ^dm,o is the present fraction of 
matter (including baryons and dark matter) and f^A.o is 
the present fraction of dark energy, we obtain 



02 
"0 



n 



A.O 



1/3 



0.677, 



(41) 



where we have used £Ib,o = 0.0455, Qdai.o = 0.1915, 
^m,o = 0.237, and fl^.o = 0.763 [l[. Introducing the 
cosmological length l\ = 4.38 1 26 m (see Appendix IB"]) . 
we obtain 



02 

/a 



"m,0 
^A,0 



1/3 / 3Q A ,, 
I 8?r 



1/2 



= 0.204. 



The value of the characteristic length 02 is 
a 2 = 8.95 10 25 m. 



(42) 



(43) 



It corresponds to the typical radius of the universe at the 
transition between the matter era and the dark energy 
era. 

Regrouping these results, we obtain the basic equations 
of the model 



P = -Pac 



PA 1 A q Pa 

~J> 9 = H 1_3 7 



Pa 



^] +1 

a 



(44) 



(45) 



(46) 



The temperature T = T„ is constant. Finally, we note 
that K = (Stt/S) 1 / 2 ^ 1 where we have introduced the 
cosmological time £a = 1-46 10 18 s (see Appendix IB"]) . 

As the universe expands from a — to a = +00, the 
density decreases from +00 to pa, the equation of state 
parameter w decreases from to —1, and the deceleration 
parameter decreases from 1/2 to —1. 

Remark: The equation of state |44|) used to describe 
the late universe dominated by dark matter and dark 
energy corresponds to a constant negative pressure (the 
velocity of sound c s = 0). This equation of state is related 



to, but different from, the linear equation of state p = 
—pc 2 . This linear equation of state describes only the 
dark energy era while the constant equation of state (|4"4"]) 
describes both the matter era and the dark energy era in 
a unified manner. Furthermore, it can be given a scalar 
field representation (see Sec. IVIII|) . It turns out that 
the equation of state (|4"4"]) is equivalent to the standard 
ACDM model (see Sec. IVIIBp . although this was not 
obvious a priori. The Chaplygin gas model @ is often 
presented as a unification of dark matter and dark energy. 
We note that a constant (not linear) negative pressure 
(|4"4"]) does the same job. In our opinion, this gives a new 
light to the ACDM model. 



B. The inflationary phase 

During the inflationary phase (a 3> 02), the density 
has an approximately constant value p ~ pa- The corre- 
sponding Hubble parameter is 



H 



'8ttG 



V 3 



-PA 



1/2 



1/2 



1 



(47) 



Numerically, p = 7.02 10" 24 g/m 3 and H = 
1.9810" 18 s _1 . Integrating Eq. (g7J), we find that the 
scale factor increases with time as 



a(t) oc e ( 87r / 3 ) 1/2t /'\ 



(48) 



This corresponds to de Sitter's solution. The radius 
of the universe increases exponentially rapidly on a 
timescale of the order of the cosmological time tA = 
1.46 10 18 s. The universe exists at any time in the future 
and there is no singularity. 



C. The matter era 

Before the late inflation (a <C 02), the universe is in the 
matter era. In that case, due to the smallness of pa, we 
can take p = and we recover the standard Einstein-de 
Sitter (EdS) model. The density is related to the scale 
factor by p ~ p\a?,/a 3 . The Fricdmann equation Q 
becomes 



H 



a ( 8irG pa«2 



1/2 



1/2 



1 /a 2 \ 3 / 2 



t A 



yielding 



- - (6-) 1/3 f 
ai V ia 



2/3 



(49) 



(50) 



where we have determined the constant of integration 
such that a = at t = 0. This is not quite correct 
because we should match the matter era to the radiation 
era (see Sec. IVTI|) . but this is a good approximation. 
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Using the approximate expression ([50]) of the scale factor 
in the matter era, we find that the time at which a = a 2 
is 



t 2 



t A 



6tt 



(51) 



Numerically, t 2 = 0.230i A = 3.37 10 17 s. This marks the 
end of the matter era and the beginning of the dark en- 
ergy era. This numerical value will be revised in Sec. IV Dl 
using the exact solution (|53p of the Fricdmann equations. 
We also have 



Pa 



(a/a 2 f 



PA 

()7T 



(52) 



During the matter era, the density decreases algebraically 
as the universe expands. 



D. The general solution 



The equation of state (|44| interpolates smoothly be- 
tween the matter era described by a density p cx a -3 
and the dark energy era described by a constant density 
p = p\ (cosmological density). It provides therefore a 
unified description of the late universe. For the equa- 
tion of state (H3J , the general solution of the Friedmann 
equation ([9]) is 



a 

<I2 



sinh 



2/3 



'67T- 



t 



The density evolves as 
P = 



PA 



tanh z ( V6tt^- 



The time at which a = a 2 is 
1 



*2 = 



argsinh(l)<A- 



(53) 



(54) 



(55) 



Numerically t 2 = 0.203t A = 2.9710 17 s. The correspond- 
ing density is p 2 = 2p A = 1.40 1CT 23 g/m 3 . For t < t 2 , 
the universe is in the matter era. Its radius increases 
algebraically according to Eq. (|50| while its density de- 
creases algebraically according to Eq. (|52"|) . For t >• t 2 , 
the universe is in the dark energy era. Its radius increases 
exponentially according to Eq. (|4"5)) while its density re- 
main approximately constant p ~ pa- From Eq. (|53[) we 
obtain the asymptotic behavior 



a(t) -W 8 ^ 1 ^ 4 - 4 ^, 



with 



1/2 



hn2 + \n( l -± 

3 \a 2 



tA- 



(56) 



(57) 
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FIG. 7. Evolution of the scale factor a with the time t in 
logarithmic scales. This figure clearly shows the transition 
between the matter era (EdS) and the dark energy era (de 
Sitter). In the dark energy era, the radius increases exponen- 
tially rapidly on a timescale of the order of the cosmological 
time tA = 1.46 10 18 s. This corresponds to a phase of late 
inflation. The universe is decelerating for a < a c and ac- 
celerating for a > a c with a c = 7.11 10 25 m. The transition 
between the matter era and the dark energy era takes place at 
d2 = 8.95 10 25 m. Coincidentally, the present universe turns 
out to be close to the transition point (eto ~ 012). 




FIG. 8. Evolution of the scale factor a with the time t in 
linear scales. 



Numerically t f = 0.708t A = 1.0310 18 s. For t > t f , 
the radius of the universe is larger than the cosmolog- 
ical length l A = 4.38 1 26 m. 

Using the results of Sec. IIV A[ we find that the universe 
is decelerating for a < a c (i.e. p > p c ) and accelerating 
for a > a c (i.e. p < p c ) where 



Pc = 3p A , a c = (l/2) 1 / 3 a 2 . 



(58) 



The time t c at which the universe starts accelerating is 
given by 



tr 



l argsinh (72^ 



(59) 
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FIG. 9. Evolution of the density p with the time t in logarith- 
mic scales. During the late inflation, the density remains ap- 
proximately constant, with the cosmological value p m in = Ph 
representing a lower bound. 




FIG. 10. Evolution of the density p with the time t in lin- 
ear scales. General relativity (cosmological constant) lim- 
its the decay of the density to the cosmological value p\ = 
7.02 10" 24 g/m 3 . 



This corresponds to the time at which the curve a(t) 
presents an inflexion point. Numerically, p c = 3 p\ = 
2.11 10~ 23 g/m 3 , a c = 0.794a2 = 7.11 10 25 m, and t c = 
0.152 t A = 2.2210 17 s. 

The evolution of the scale factor and density as a func- 
tion of time are represented in Figs. [71 1101 in logarithmic 
and linear scales. 



E. The present universe 

The Hubble radius, the density and the Hubble time 
of the present universe are a = 0.302?a = 1.32 10 26 m, 
Po 



1.31/5A = 9.20 1 0- 24 g/m 3 , and # ~ 



0.302t A 



4.41 10 17 s. Using Eq. (|5"3"|) . the age of the universe is 



given by 



to 



argsinh 



3/2' 



(60) 



Numerically, t = 0.310£a = 4.5410 17 s. According to 
Eq. (j45|) , the present values of the deceleration parameter 
and of the equation of state parameter are 



l-3fi 



A.O 



w 



-^A,0- 



(61) 



Numerically, <jo = —0.645 and = —0.763. 

It is striking to note that the present size of the uni- 
verse ao = 1.32 10 26 m is of the order of the scale a-i — 
8.95 1 25 m (a = 1.48a 2 ). Therefore, we live just at the 
transition between the matter era and the dark energy 
era (see bullets in Figs. I71ITU]). Another way to state this 
result is to say that the present ratio f2A,o/^?n,o = 3.22 
between dark energy and dark matter is of order unity. 
This coincidence is intriguing and often referred to as the 
"cosmic coincidence problem" [Io| . Several theories have 
been proposed to explain why Q\fi/fl m _Q ~ 1 fill ]. How- 
ever, this may be just a pure coincidence without deeper 
reason. Life (and researchers inquiring about cosmol- 
ogy) may have emerged ~ 14 Gyrs after the Big Bang, 
precisely at the epoch where fiA/f2 m ~ 1. We leave this 
"problem" open. 



A MODEL OF PECULIAR CYCLIC 
UNIVERSE 



VI. 



The equation of state (fT0|) with n < and k > 
leads to a model of non-inflationary universe exhibiting 
a future peculiarity. We take the same values of a, n, p*, 
and a,2 as in the previous section. We also take k = p\. 
The basic equations of the peculiar model are 



w 



P = Pac , 



Pa 1 , 

7' q= 2 { 



3 PA 



P = PA 



- 1 



(62) 



(63) 



(64) 



The temperature T = T* is constant. As the universe 
expands from a — to 02, the density decreases from 
+00 to 0, the parameter w increases from to +00 and 
the deceleration parameter increases from 1/2 to +00. 

For the equation of state (f6"2")h the general solution of 
the Friedmann equation (|9]) is 



a 

(12 



sm 



2/3 



t 



'6ir — 



t A 



(65) 



11 




FIG. 11. Evolution of the scale factor a with the time t. The 
dotted line corresponds to the real evolution of the universe 
according to the model of Sec. [V] The bullet corresponds 
to the value of the scale factor at the time to = 0.310tA = 
4.54 10 17 s corresponding to the age of our universe. At this 
epoch, the cyclic universe has deviated relatively strongly 
from the real evolution. 
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a, 
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FIG. 12. Evolution of the density p with the time t. The den- 
sity vanishes periodically at the moment where the universe 
reaches its maximum size. On the other hand, it becomes 
infinite when the universe contracts to a point (Big Bang or 
Big Crunch). 



The density evolves as 

P_ = 
PA 



tan' 



(66) 



This solution describes a peculiar cyclic universe. The 
universe reaches its maximum radius a 2 = 0.204Za = 
8.95 10 25 m for the first time at 



ti = — 
24 



1/2 



(67) 



Numerically, t 2 = 0.362t A = 5.2910 17 s. For t < t 2 , we 
can make the approximation p ~ 0. Therefore, this uni- 
verse evolves just like the Einstein-de Sitter (EdS) uni- 
verse describing the matter era (see Sec. IV C|) . Its radius 



increases algebraically according to Eq. (|50|) while its 
density decreases algebraically according to Eq. (|52"|) . 
When t — > t 2 , the radius approaches its maximum value 
according to 



a 2n , 

i--~ w te-ty 



Correspondingly, the density decreases to zero as 
P 67r 



PA 



-(h-ty. 



(68) 



(69) 



For t > t 2 , the radius decreases with time while the den- 
sity increases (this corresponds to the solution of Eq. 
(fT4"|) with e = —1). The radius vanishes at t = 2t 2 
(Big Crunch) , then increases again. The density vanishes 
when a = a 2 and is infinite when a = 0. This process 
continues periodically. Although there is no divergence 
at t = t 2 , this model is peculiar because, when the ra- 
dius reaches its maximum value a 2 , the density vanishes 
(p = 0), so the universe becomes empty. In a sense, the 
universe "disappears" at t = t 2l before re-appearing at 
t > t 2 , and this in a periodic manner. According to this 
model, the universe would disappear for the first time 
at t 2 = 0.362iA = 5.29 10 17 s that is not very far from 
the present age of the universe to = 0.310£a = 4.54 10 17 s. 
Therefore, if we were living in this cyclic universe, the end 
of the world (in the sense of the whole universe) would 
take place in about 2.38 billion years. We could fear to 
live in this universe. However, this universe is always de- 
celerating in the phase of expansion. Therefore, the fact 
that our universe is accelerating (an observation that was 
made only recently @) shows that we do not live in this 
cyclic universe. In a sense, this is reassuring. 

The evolution of the scale factor and density as a func- 
tion of time are represented in Figs. HT1and[T2l 

Remark: The equation of state (|44|) corresponds to a 
constant positive pressure (the velocity of sound c s = 
0). This model turns out to be equivalent to the cold 
dark matter model with a negative cosmological constant 
(anti- ACDM model) , although this was not obvious a pri- 
ori. This is an example of cyclic universe corresponding 
to k = 0, A < and p = 0. Other types of cyclic uni- 
verses (corresponding to k = ±1, A < and p — 0, or 
k = +1, A = and p = 0) were studied by Friedmann 
[t2T |. Einstein [t|, and Lemaitre [3, and were called 
"phoenix-universes" by Lemaitre. 



VII. A SIMPLE MODEL FOR THE WHOLE 
EVOLUTION OF THE UNIVERSE 

A. The early universe 

In Paper I, we have described the evolution of the 
early universe, corresponding to the transition between 
the pre-radiation era and the radiation era, by a single 
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equation of state of the form (flOl) with a = 1/3, n = 1 
and k = — 4/(3pp), namely 



1 



P = oP(! - ±p/PpV 



(70) 



This equation of state leads to a non-singular early in- 
flationary universe. The density is related to the scale 
factor by 

PP 



P = 



(71) 



(0/0!)* + 1' 

where a\ = 2.61 10 _6 m. This characteristic scale marks 
the transition between the pre-radiation era and the ra- 
diation era. 

When a<(ii, the density tends to a maximum value 

P = Pmax=PP, (72) 

identified with the Planck density pp = 5.16 1 99 g/m 3 . 
This leads to a phase of early inflation 

aaeW 3 ' 1 ' 2 *, (73) 

whose timescale is the Planck time tp = 5.39 10 _44 s. 

When d>ai,we recover the equation p ra d ~ ppaf/a 4 
corresponding to the pure radiation described by an 
equation of state p = pc 2 /3. The conservation of p ra .a,a A 
implies that p P a\ = p r ad,oa^. Writing p ra d,o = &rad,oPo, 
we can rewrite Eq. (|7ip in the form 



ClradfiPO 



(74) 



(a/ao) 4 + (ai/a ) 4 ' 
where a±/ao = 1.97 1 - 32 . When a ^> a±, it reduces to 

Clrad, OpO 



Prad 



(a/ao) 4 



(75) 



which corresponds to the pure radiation. In the radiation 
era, the scale factor, the density, and the temperature 
evolve as a oc t 1 / 2 , p oc t~ 2 and T oc i" 1 / 2 . 

The Friedmann equation ^ corresponding to Eq. ([74]) 
can be written 



H 



(a/a ) 4 + (fli/ao) 4 



(76) 



It has the analytical solution given by Eq. (1-128). The 
transition between the pre-radiation era and the radi- 
ation era takes place at a typical time t\ = 23.3tp = 
1.25 10~ 42 s corresponding to a = a\. This time also cor- 
responds to the inflexion point of the curve a(t). The 
universe is accelerating for t < t\ and decelerating for 

t > h. 

Remark: The equation of state (|70p smoothly interpo- 
lates between a negative pressure vacuum energy regime 
(p = —pc 2 , p = pp) and a relativistic radiation phase 
(p = pc 2 /3, p oc a~ 4 ). This amounts to summing the 
inverse of the densities of these two components taken 
independently. Indeed, Eq. (|74fl can be rewritten as 



1 

P 



1 



Prad 



l 

pp' 



(77) 



B. The late universe 

In this paper, we have described the evolution of the 
late universe, corresponding to the transition between the 
matter era and the dark energy era, by a single equation 
of state of the form (|T0|) with a = 0, n — — 1 and k — 
— pa, namely 



P 



-PAC 



(78) 



This equation of state leads to a non-singular late in- 
flationary universe. The density is related to the scale 
factor by 



P = PA 



an 



1 



(79) 



where a 2 = 0.204Za = 8.95 1 25 m. This characteristic 
scale marks the transition between the matter era and 
the dark energy era. 

When o> a2, the density tends to a minimum value 



P 



PA, 



(80) 



identified with the cosmological density p\ = 
7.02 10~ 24 g/m 3 . This leads to a phase of late in- 
flation 



3 (87T/3) 1 / 2 t/tA 



(81) 



whose timescale is the cosmological time t\ = 1.46 10 18 s. 

When a <C 02, we recover the equation p m ~ pAa^/a 3 
corresponding to the pressureless matter described by an 
equation of state p = 0. The conservation of p m a 3 implies 
that p A »2 = Pm.oao- Writing p m . = fi m , Po and p Ai0 = 
^A.oPOi we can rewrite Eq. (|79p in the form 

Qm,oPo 



(82) 



(83) 



(a/ao) 3 
When a <C a 2 , it reduces to 

Pm — 7 / r; 
(a/a J 3 

which corresponds to the pure matter. In the matter era, 
the scale factor and the density evolve as a cx i 2 / 3 and 
p (X t~ 2 . 

The Friedmann equation ((9]) corresponding to Eq. (|82|) 
can be written as 



H 
Ho 

or, equivalently, as 



TT = \/fim,o [( fl o/a) 3 + (a /a2) 3 ] 

tin V 



H 
Ho 



CI 



A,0) 



(84) 



(85) 



with fl m + CI \ = 1. It has the analytical solution given 
by Eq. (|53l) . This can be rewritten as 

1/3 

sinh 2/3 f ^^CJ^Hot) . (86) 



a 



- - / ; 



CI 



A.I) 
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The evolution of the density is given by 



baryonic matter 



P_ 

Po 



^A,0 



tanh 2 (fvT^o-ffo*) 



(87) 



The transition between the matter era and the dark en- 
ergy era takes place at a typical time t% = 0.203tA = 
2.9710 17 s at which a — 02. The universe is decel- 
erating for t < t c and accelerating for t > t c where 
t c = 0.152t A = 2.2210 17 s. Using the results of Sec. 
and Appendix [B] we can write 



02 

ao 



^A,0 



1/3 



ao 



( "m,0 

\2r2A,o 



1/3 



1 2 1 



iargsinh(l), 



A.O 



1 2 1 



#0 3 y/fl A ,o 



argsm 



inh(l/V2). 



(89) 



(90) 



Numerically, 02 = 0.677ao and a c = 0.538ao. The age of 
the universe is given by Eq. (|60|) which can be rewritten 



1 2 
to = ~rr~ 77 



1 



Ha 3 y/^a.o 



argsinh 



m,0 



1/2' 



(91) 



Numerically, t = 0.310i A = 4.53 10 i7 s. 

Remark: The equation of state (|78[) smoothly inter- 
polates between a non-relativistic matter phase {p = 0, 
p oc a~ 3 ) and a negative pressure dark energy regime 
(p = —pc 2 , p = p\). This amounts to summing the 
density of these two components taken independently. 
Indeed, Eq. (|52")l can be rewritten as 



P = Pn 



PA- 



(92) 



Therefore, the equation of state (fT8"f is equivalent to the 
standard ACDM model [l|. 



C. Connection between the early and the late 
universe 

The previous sections reveal a deep "symmetry" be- 
tween the early and the late evolution of the universe. In 
the early universe (n = 1), we have to add the inverse of 
the densities, and in the late universe (n = —1) we have 
to add the densities themselves. This suggests that the 
universe must be described by two equations of state. 

In order to obtain a single equation describing the 
whole evolution of the universe, we just have to add the 
contribution of radiation (including pre-radiation) 



Prad 



Qrad,0PQ 



(a/a ) 4 + (ai/ao) 4 



(93) 



dark matter 



and dark energy 



Pb = 



Pdm 



^BfiPO 

(a/a ) 3 ' 

^DMfiPO 

(a/a ) 3 



PA = ^A,0P0> 



(94) 



(95) 



(96) 



in the Friedmann equation (|TJ]). The radiation dominates 
in the early universe, while baryonic matter, dark matter 
and dark energy dominate in the late universe. Now that 
the radiation phase has been regularized in order to avoid 
the initial singularity, it is necessary to state explicitly 
that matter is formed after radiation^. To that purpose, 
it is sufficient to add a Hcaviside function H(a — 01) in 
the matter density. If we write Q* B = Qb,oH(ci — a\) 
^DM.oH(a — ai), the Friedmann equation 



and ft 



DM,0 



© can be written 



H 
#0 



k/.O 



"s,o 



n 



(a/a ) 4 + (ai/a ) 4 (a/a ) 



+ 



DM,0 



^A,0, 



(a/a ) 3 

(97) 

with Qb + ^rad + &dm + = 1- For cl\ = 0, we re- 
cover the standard model of the universe [![. It exhibits 
a singularity at t — (Big Bang). For a\ ^ 0, we obtain 
a regularized model without primordial singularity. The 
universe always existed in the past but, for t < 0, it has 
a very small radius, smaller than the Planck length. At 
t = 0, it undergoes an inflationary expansion in a very 
short lapse of time of the order of the Planck time and 
connects to the standard model. A nice feature of this 
model is its simplicity since it incorporates an "inflation 
phase" in a very simple and very natural manner. We just 
have to add a term +{a\/ao) A in the standard equation 
of H/Ho given in [l[. Therefore, the modification implied 
by Eq. (|97[) to erase the initial singularity is very natu- 
ral. On the other hand, the equation of state (|75)l in the 
late universe is equivalent to adding the contribution of 
matter and dark energy individually, as in the standard 
model [l|. This does not bring any modification to the 
usual equation, which is a virtue of this description since 
the standard model works well at late times. 



1 In the standard model, the radiation density diverges as a -4 
when a — > while the matter density diverges as a -3 which is 
subdominant. Therefore, we can safely extrapolate the matter 
density to a = (its contribution is negligible anyway). In the 
present model, the radiation density J93P does not diverge any- 
more. Therefore, it becomes necessary to state explicitly that 
matter is formed after radiation in order to avoid spurious diver- 
gences of the matter density when a — ¥ 0. 
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FIG. 13. Temporal evolution of the scale factor in logarithmic 
scales. The universe exists at all times in the past and in the 
future and there is no singularity (aioniotic universe). The 
early universe undergoes a phase of inflation that takes it to 
the radiation era. This is followed by the matter era and the 
dark energy era responsible for the accelerated expansion of 
the universe. The universe exhibits two types of inflation: An 
early inflation corresponding to the Planck density pp due to 
quantum mechanics and a late inflation corresponding to the 
dark energy density pA due to the cosmological constant (gen- 
eral relativity). The evolution of the early and late universe 
is remarkably symmetric. We have represented in dashed line 
the standard model leading to a primordial singularity (Big 
Bang). The dotted lines corresponds to the model of Sec. 
IVIIBI where the radiation is neglected. This analytical model 
provides a good description of the late universe. We have also 
represented the location of the present universe that is just 
at the transition between the matter era and the dark energy 
era. 



FIG. 15. Temporal evolution of the temperature of the radia- 
tion in logarithmic scales. The temperature increases consid- 
erably during the early inflation and decreases considerably 
during the late inflation. The present value of the tempera- 
ture is ~ 2.7K. 



yields 



a/a 



dx 



Hot, 



x 4 +(ai/a ) 4 



11 DM,0 
-=3 



o 



A.O 



(98) 

where at = ai = in the standard model and <Xj = lp in 
the regularized model. The age of the universe is 



1 



dx 



^DM,0 
— 3— 



(99) 
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FIG. 14. Temporal evolution of the density in logarithmic 
scales. The density goes from a maximum value p max = pp 
determined by the Planck constant (quantum mechanics) to 
a minimum value p m i„ = pa determined by the cosmological 
constant (general relativity) . These two bounds are responsi- 
ble for the early and late inflation of the universe. In between, 
the density decreases as t~ 2 . 



The evolution of the scale factor with time is obtained 
by solving the first order differential equation (|97p . This 



Of course, for the determination of the age of the uni- 
verse, we can neglect the pre-radiation era and take 
ai = (strictly speaking, the age of the universe is infi- 
nite since it has no origin; however, we define the age of 
the universe from the time t = at which a = lp). We 
obtain the standard result to = 1-03 H^ 1 = 4.53 10 17 s = 
14.4 GyrQ. Actually, we find the same result if we neglect 
radiation and use the analytical expression (|91[) instead 
of Eq. (|99|) . It is rather fortunate that the age of the 
universe almost coincides with the Hubble time Hq 1 . 

In Figs. [T3TiT5l we have represented the evolution of 
the radius, density and temperature of the universe as 
a function of time. The universe exhibits two types of 
inflations: An early inflation due to the Planck density 
Pp = 5.16 10 99 g/m 3 and a late inflation due to the cosmo- 
logical (dark energy) density p\ = 7.02 10 _24 g/m 3 . The 
early inflation can be described by a polytropic equation 
of state with positive index n = 1. The late inflation can 



2 The Hubble constant is usually written as Ho = 
2.268 fi7 10 — 18 s _1 where the dimensionless parameter /17 is 
about 10% of unity [l|. For simplicity, we have taken hj = 1 in 
the numerical applications. The current value is /17 = 1.05±0.05. 
If we take hi = 1.05, the age of the universe is to = 13.7 Gyr. 
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be described by a polytropic equation of state with neg- 
ative index n = — 1. There exists a striking "symmetry" 
between the early and the late evolution of the universe, 
the cosmological constant in the late universe playing the 
same role as the Planck constant in the early universe. 
In particular, Fig. [14] shows that the density varies be- 
tween two bounds p ma x = Pp and p m i n = p\ that are 
fixed by fundamental constants (see Appendix [B]). These 
values differ by a factor of the order 10 122 . The early 
universe is governed by quantum mechanics (h) and the 
late universe by general relativity (A). 
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FIG. 16. Evolution of the scale factor in the singular model 
corresponding to k < 0. The universe starts at t = from a 
singularity at which the radius is finite but the density and 
the temperature are infinite. This is followed by the radiation 
era, the matter era, and the dark energy era. 



For the sake of completeness, we have also represented 
in Fig. [in] the evolution of the scale factor in the model 
with the new type of initial singularity considered in Pa- 
per I. This model is obtained by replacing the sign + 
by — in Eq. (|93[) and in the following formulae (we 
must then take a% = a\ in Eq. ([95]) ). In that case, 
the universe emerges "from nothing" at t = 0. It has 
an infinite density and temperature but a finite radius 
ai = 2.61 10 _6 m. This universe becomes physical only 
for t > U = 0.154t P = 8.32 10 _45 s when the velocity of 
sound becomes smaller than the speed of light. 

Remark: An additional refinement can be introduced if 
dark matter is made of BECs instead of being pressure- 
less. In that case, the term fl* DM /(a/ao) 3 should be 
replaced by Q* DM . /[(a/a ) 3 T (aBEc/ao) 3 } as proposed 
in d, Qj} . It is shown in these papers that the presence of 
BECs can substantially accelerate the formation of the 
large scale structures of the universe. 



D. The most natural universe 

We may wonder whether the actual evolution of the 
universe can be predicted, without making any observa- 
tion, from the principle of "simplicity" and "aesthetics" 
that was defended by great scientists such as Einstein, 



Eddington, and Chandrasekhar [18j. We make the fur- 
ther assumption that the universe is non-singular and 
non-peculiar. Therefore, we try to construct the simplest 
non-singular model of universe. 

The Friedmann equations ©-([1]) governing the evolu- 
tion of the universe require the value of the curvature of 
space fc, the cosmological constant A, and the equation 
of state p = p(p). Concerning the curvature of space, the 
most natural guess is k = (this leads to the simplest 
mathematical models). Actually, this value is consistent 
with observations. Even if an apparently flat universe 
could be the result of the primordial inflation (leading to 
k/a 2 <C 1) 0: h ma y weu be that k = exactly. This is 
what we assume. The cosmological constant can be elim- 
inated from the basic equation since it can be taken into 
account in the equation of state. Concerning the equation 
of state, the simplest guess is to take p = 0. This choice 
was made by early cosmologists. A pressureless universe 
describes the matter era (EdS model). However, this 
equation of state neglects relativistic effects which be- 
come important in the past when the universe was hot. 
Therefore, cosmologists considered a linear equation of 
state of the form p = ape 2 . For a = 1/3, this equation of 
state describes the (relativistic) radiation era. However, 
this equation of state is not fully satisfactory since it leads 
to a primordial singularity (Big Bang). The early infla- 
tion and the late inflation (dark energy) are usually de- 
scribed by a linear equation of state of the form p = —pc 2 
leading to a constant density. However, this equation of 
state does not contain in itself the value of the constant 
density (a feature that we regard as a drawback). On 
the other hand, being linear, it cannot be "added" to an- 
other linear equation of state p = ape 2 to describe the 
transition between the pre-radiation era (a = —1) and 
the radiation era (a = 1/3) in the early universe, or the 
transition between the matter era (a = 0) and the dark 
energy era (a = — 1) in the late universe. 

The next simplest choice is therefore to consider a 
mixed equation of state of the form p = (ap + kp 1+1 / n )c 2 
with a linear component and a polytropic component. It 
comes naturally from the polytropic model that the in- 
dices n > describe the early universe while the indices 
n < describe the late universe. This suggests that the 
early and the late universe must be described separately 
by two different equations of state. We must take a = 1/3 
(radiation) in the equation of state describing the early 
universe (n > 0) and a = (matter) in the equation of 
state describing the late universe (n < 0). By analyzing 
all the possibilities, we have found that only the models 
with k < can produce a non-singular (or non-peculiar) 
universe. The models with k > lead to past or fu- 
ture singularities/peculiarities. Therefore, the polytropic 
pressure must be negative. Then, it comes naturally from 
the polytropic model that the density has a maximum 
value Pmax in the early universe (for n > 0) and a min- 
imum value Pmin in the late universe (for n < 0). It is 
natural to identify the upper bound with the Planck den- 
sity pp (quantum mechanics) and the lower bound with 
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the cosmological density p\ (general relativity) . This de- 
termines the constant k which is k = — 4/(3/3p™) in the 

early universe and k = — 1/ Pa in the late universe. Fi- 
nally, according to the principle of simplicity, it is natural 
(but not quite compulsory) to select the index n = +1 
(quadratic equation of state) in the early universe and 
the index n = — 1 (constant equation of state) in the late 
universe. This leads to the model represented in Figs. 
1 1 341 1 51 featuring two symmetric phases of inflation. This 
evolution could have been predicted without making any 
observation of the universe, by only considering the sim- 
plest equation of state that does not yield any singularity 
or peculiarity. Interestingly, the observations of the real 
universe are remarkably consistent with that model. 

In conclusion, the most natural non-singular model of 
universe corresponds to a solution of the Einstein equa- 
tions that exists eternally in the past and in the future. 
This has been called aioniotic universe in Q since it 
has no beginning nor end. It is the only model of the 
family of equations of state pU|) that has no singular- 
ity/peculiarity. For that reason, it may be selected by 
Nature. Furthermore, this universe is very "symmetric" 
in its early and late phased). This symmetry can be seen 
in the values of the polytropic indices n = +1 and n = — 1 
that characterize the early and the late universe, respec- 
tively. All the formulae obtained in Paper I and in the 
present paper are symmetric with respect to each other 
(compare Sees. IVII Al and IVII B|) . This symmetry is also 
apparent in the two types of inflation that describe the 
early and late universe: An early inflation due to the 
Planck density pp (quantum mechanics) and a late in- 
flation due to the cosmological density p\ (general rela- 
tivity). They correspond to the upper and lower bounds 
of the density shown in Fig. Q3J which is also strikingly 
symmetric. Therefore, this universe is the simplest and 
most symmetric non-singular cosmological solution of the 
Einstein equations. 

Of course, some reservations should be made: The uni- 
verse may not be "simple" . In particular, the very early 
universe (before the Planck time) may not be described in 
terms of an equation of state p{p), or even in terms of the 
Einstein equations, as we have assumed. The description 
of the very early universe may require the development 
of a theory of quantum gravity that does not exist for 
the moment. An interesting description of the early in- 
flation has been given by [19j in terms of a quantized 
model based on a simplified Wheeler-DcWitt equation. 
In that model, a quantum tunneling process explains the 
birth of the universe with a well defined size after tunnel- 
ing. Therefore, other inflationary scenarios are possible 
in addition to the one given in Paper I in terms of the 



3 We may argue that a model of late universe with pa = (cor- 
responding to a vanishing cosmological constant) is also non- 
singular and simpler. However, if we add a principle of "symme- 
try" (~ aesthetics) . in addition to the principle of "simplicity" . 
we are left with the model discussed previously. 



generalized equation of state (|10j) . However, our aim was 
to explore all the consequences of this generalized equa- 
tion of state, pushing it to its limits. 



E. A cosmological constant "problem"? 

The cosmological constant A is equivalent to a con- 
stant density pa = A/(8irG) called dark energy. Its 
value resulting from observations is p\ = 7.02 10 _24 g/m 3 
(see Appendix IBjl . Since the late universe is very di- 
lute, it is oftentimes argued that the cosmological den- 
sity pa should correspond to the vacuum energy den- 
sity due to quantum fluctuations. However, according 
to particle physics and quantum field theory, the vac- 
uum energy density is of the order of the Planck density 
Pp = 5.16 10"g/m 3 which is 10 122 times larger than the 
cosmological density. This leads to the so-called cosmo- 
logical constant problem [2(| . 

Actually, as illustrated in Fig. HH the Planck density 
and the cosmological density represent fundamental up- 
per and lower density bounds acting in the early and late 
universe, respectively. It is not surprising therefore that 
they are so different: pk^ Pp- Because of these bounds, 
the universe undergoes two phases of inflation. The infla- 
tion in the early universe is due to quantum mechanics 
and is related to the Planck density pp. The inflation 
in the late universe is due to the cosmological constant 
and is related to the cosmological density p\. Quantum 
mechanics is negligible in the late universe. The cosmo- 
logical constant should be interpreted as a new funda- 
mental constant of physics. It applies to the very large 
universe (cosmophysics) exactly like the Planck constant 
applies to the very small universe (microphysics). Ac- 
tually, there is a complete symmetry between the small 
and large universe where Ti and A play symmetric roles. 
Therefore, we interpret the cosmological constant as a 
fundamental constant of physics describing the cosmo- 
physics (late universe) in the same sense that the Planck 
constant describes the microphysics (early universe). 

If this interpretation is correct, the origin of the dark 
energy density p\ should not be sought in quantum me- 
chanics, but in pure general relativity. In this sense, the 
cosmological constant "problem" may be a false problem. 
If A is a fundamental constant of physics, independent 
from the others, its value should not cause problem. It 
is fixed by nature, just like the value of G, c, and h. One 
can just expect that pp is "very large" and p\ is "very 
small". Of course, the origin of the cosmological con- 
stant still needs to be understood by developing a theory 
of cosmophysics. In addition, it is extremely important 
to understand why pp and pa represent upper and lower 
bounds, and if these bounds are as fundamental as, for 
example, the bound on the velocity fixed by the speed of 
light. 

Remark: We may note that the appearance of maxi- 
mum and minimum density bounds from the polytropic 
equation of state pop is, in some sense, similar to the 
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appearance of the maximum mass of relativistic white 
dwarf stars discovered by Chandrasekhar 21[ also from 
a polytropic equation of state. 



VIII. SCALAR FIELD MODELS 



The phase of inflation in the very early universe is usu- 
ally described by a scalar field |3|. On the other hand, 
in alternative theories to the cosmological constant, the 
present-day acceleration of the universe (dark energy) 
is also described by a scalar field called quintessence 
@, 0. A tachyon field @, [l| has also been considered 
due to its connection with string theory [l6| . Actually, 
to any "fluid" described by a barotropic equation of state 
p = p(p), it is possible to attach a scalar field with a po- 
tential V(<f>). In this section, we determine the potential 
of the scalar field (quintessence and tachyon field) corre- 
sponding to the equation of state (fTU|) using the general 
methodology exposed in [(|. We treat polytropic indices 
n and polytropic constants k of arbitrary sign but we 
assume a + 1 + fcp 1 /™ > (i.e. w > — 1). The case 
a + 1 + kp x l n < (i.e. w < —1) leads to phantom 
(ghost) fields that are studied in Paper III. 



A. Quintessence 

Quintessence [H, 0] is described by an ordinary scalar 
field <p minimally coupled to gravity. The scalar field 
evolves according to the equation 



form 



3H<p+ — = 0, 



(100) 



where V((f>) is the potential of the scalar field. The 
quintessence scalar field tends to run down the potential 
towards lower energies. The density and the pressure of 
the universe are related to the scalar field by 

pc 2 = \ft + V(<f>), P = ^4?-v{4>). (101) 

When the kinetic energy dominates the potential energy, 
we obtain the equation of state p — pc 2 of stiff matter. 
When the potential energy dominates the kinetic energy, 
we obtain the vacuum equation of state p = —pc 2 . 
From Eq. ([TOT]) , we get 



4> 2 = (1 + w)pc 2 



(102) 



where we have written p = wpc 
in order to have 4> 2 > 0. Using <f> 
Fricdmann equation (|9]) valid for a flat universe, we get 



We assume w > — 1 
(d<j)/da)Ha, and the 



da 



3c 2 
SttG 



1/2 



(103) 



For the equation of state (fTU|). using Eqs. (fTTj) and (1-69), 
and setting R = a/a*, we can rewrite Eq. (|103[) in the 



d<t> _ f 3e 2 V /2 y^TT fl3(i+°)/2" 



With the change of variables 



fl 3(a+l)/2n ^ = (, 



3^2" 



2n 



we find that 



dx 



0, (105) 



(106) 



Vx 2 =F 1 
For k > 0, we obtain 

tp = Argcosh(x), (107) 
and for k < 0, we get 

%j) = Argsinh(a;). (108) 
On the other hand, according to Eq. (|10ip . we have 
1 



V = ^(l-w)pc 2 . 



(109) 



For the equation of state (JTUJ), using Eqs. (fTTj) and (1-69), 
we obtain 



V = -p*c 



1_ 2 (l-a)x 2 T2 



2^ (x 2 =fl) n+1 



(110) 



For k > 0, since x = cosh(-0) the scalar field potential is 
explicitly given by 



1 2 (l-a)sinh 2 ?A-(a + l) 
V{tp) = t;P*c ~ 



sinh 2( ' Tl+1 '' -0 



(111) 



For k < 0, since x = smh(?/>) the scalar field potential is 
explicitly given by 

In these models, ij> > 0. Let us consider particular cases. 

(i) For n = 1, k = —4/ (3pp) and a = 1/3 (non-singular 
early universe), we obtain 



1 2 cosh 2 ?A + 2 

V W) = oP pC 



3 rJ cosh 4 ip 



(113) 



and i? 2 = sinh(^). 

(ii) For n = 1, k = 4/(3pp) and a = 1/3 (singular 
early universe), we obtain 



1 2 sinh ^ - 2 

= oPf c — ■ r 4 , ■ 
o smn ■(/> 



(114) 



and i? 2 = cosh(-0). 
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(iii) For n = — 1, k = — p\ and a — (late universe, 
ACDM model), we obtain 



V{ip) = -p A c 2 (cosh 2 ^ + l) 



and i?~ 3 / 2 = sinh(V')- 

(iv) For n = — 1, fc = p\ and a 
anti-ACDM model), we obtain 



VW = -p A c 2 (sinli 2 V-l) 



(115) 
(late universe, 
(116) 



and i?~ 3 / 2 = cosh(V'). 

(v) For n = —1/2, fc < and a = (Chaplygin gas), 
we obtain 

^)^c 2 (cosh^ + ^), (117) 

and i?~ 3 = sinli(^). This returns the result of Q. 

(vi) For n = — 1/2, k > and a = (anti-Chaplygin 
gas), we obtain 



Vty) 



- 2 P,c 



sinIn/> 



sinh-0 / ' 



(118) 



and i?~ 3 = cosh( , 0). 

The case a = — 1 and fc > must be treated specifi- 
cally (see Appendix A of Paper I and Appendix lAl of the 
present paper). Repeating the preceding procedure, we 
find that the potential of the scalar field is 



V{</>) = p*c 2 



1 



2 2 



\n\& 

2^G J V 1 ~ 12ttG0 2 
2nG l0 



InR = 



nc 



-, (119) 



(120) 



Finally, for a linear equation of state p = ape 2 with 
a > —1, writing the relation between the density and 
the scale factor as p/ p* = (a»/a) 3 ( 1+a ', we obtain [(|: 

V{<f>) = \p*c 2 {l - a)e- 3 ^(^) 1/2 *, (121) 



3c 1 



8ttG/ 



1/2 



\/l + a lni?, 



(122) 



where R = a/a*. Since R oc t 2 ^ 3t - 1+a ^, the scalar field 
evolves with time as <j> = [c 2 /6ttG(1 + a)] 1 / 2 Int. We 
note that V((f>) = for a — 1 (stiff matter). On the 
other hand, there is no scalar field corresponding to the 
vacuum equation of state p = ~pc 2 (i.e. a = — 1) used 
to describe inflation. By contrast, the equation of state 
(fT0"| admits a scalar field representation in all cases, in 
particular for the models with k < exhibiting a phase 
of inflation. This may be an advantage of the polytropic 
equation of state (|10[) with respect to the linear equation 
of state p = —pc 2 . 



B. Tachyon field 

A tachyon field @, [l5[ has an equation of state p = 
wpc 2 with — 1 < w < 0. This scalar field evolves accord- 
ing to the equation 



The density and the pressure are given by 



pc 



P =-V(0)dl- ( p 2 



(123) 



(124) 



From these equations, we obtain 

4> 2 = l + w, 

where we have written p = wpc 2 . Using <; 
and the Friedmann equation ©, we get 

2 \ 1/2 



_ / 3c 2 y /z yr 



da \8nGJ ^fpc^a 



(125) 
(d<j)/da)Ha, 

(126) 



For the equation of state (|TU|). using Eqs. (fTTj) and (1-69), 
we can rewrite Eq. f|126[) in the form 



d<j> 



3c 2 ^ 1/2 



V a + 1 p3(l+q)/2n 



dfl \8nG J R 

x r J R 3 ( i +«)/« T i 

With the change of variables 



(n-l)/2 



• (127) 



(128) 

we find that 



0= / (x 2 T l) ( "~ 1)/2 ^. 



(129) 



On the other hand, from Eq. (|124[) , we have 

V 2 = -wp 2 c 4 . (130) 



For the equation of state (JTUJ), using Eqs. (fTTj) and (1-69), 
we obtain 



2„4 



v = -pic 



ax 2 ± 1 
(x 2 T l) 2n+1 ' 



(131) 



Therefore, the scalar field potential V (i/j) is given in para- 
metric form by Eqs. (|129j) and (|131|) . Let us consider 
particular cases. 

(i) For n = 1, we find that x = ip. Therefore, we obtain 



V 2 



2 4 Q-0 2 ± 1 



(132) 



19 



and i? 3 ( Q+1 )/ 2 = ip. In the case k < with a < 0, where 
w is always between —1 and 0, we obtain 



V(iP) = 



\J\ — aip 2 

(^2 + 1 )3/2' 



(133) 



with ip > 0. We note that for ?i = 1, fc = — 4/(3pp) 
and a = 1/3 (early universe), the parameter w becomes 
positive when R > R w = 3 1 / 4 so that the tachyon field 
is not defined for all times. Similarly, in the case k > 0, 
since it' — > +oo when R — > 1, the tachyon field is not 
defined for all times. 

(ii) For n = — 1 and k > 0, we find that x = 
— 1/ tanh?/;. Therefore, we obtain 



V 2 = -P 2 c± 



tanh ip 



1 



sinh ip 



(134) 



and i? 3 ( a+1 )/ 2 = -tanh-0, with ip < 0. Since w +oo 
when i? — » 1, the tachyon field is not defined for all times. 
In particular, for n = — 1, k = p\ and a = (anti-ACDM 
model), the parameter w is always positive so that the 
tachyon field is not defined. 

(iii) For n = — 1 and k < 0, we find that x = tan^. 
Therefore, we obtain 



V 2 = -^ c 4 (atan 2 ^- 1) 



1 



cos 2 ip ' 



(135) 



and i?3(a+i)/2 = i/t a nV> with < ip < n/2. If a < 0, w 
remains between —1 and 0, so the tachyon field is well- 
defined. In particular, for n = — 1, k = — p\ and a = 
(late universe, ACDM model), we find 



and i? 3 / 4 = VT 
(vi) For n = - 
we find that V{(f)) 



- V> 2 /^ with < ip < 1. 
1/2, a = 0, and fc < (Chaplygin gas), 
p^c 2 is constant. The tachyon field is 



not defined in the anti-Chaplygin gas (n = —1/2, a = 0, 
and k > 0). 

The case a = — 1 and /c > must be treated specifi- 
cally (see Appendix A of Paper I and Appendix [5] of the 
present paper). Repeating the preceding procedure, we 
find that the potential of the scalar field is 



v^ 2 = P y 



xil- 



In R = sgn(n) 



27rGp*(n + l) 2 
2nGp^n + lf 



2nGp*{n+lY 



2n/(n+l) 



1 



04n/(n 



+ 1 



1/(71+1) 



02/(«+l) 



l/(n+l) 



/,2/(n+l) 



(140) 



(141) 



For n = 1 and t > i,-, and for n = — 1 and t < tf, 
corresponding to the period over which the velocity of 
sound is smaller than the speed of light, the parameter w 
is between — 1 and so that the tachyon field is perfectly 
well-defined. 

Finally, for a linear equation of state p = ape 2 with 
— 1 < a < 0, writing the relation between the density 
and the scale factor as p/p* = (a*/a) 3 ( 1+a \ we obtain 



v{4>) = f 



6ttG, 



(142) 



PAC 

COS'0' 



(136) 



and R?f 2 = 1/tan^ with < V < tt/2. 

(iv) For n = —2 and k > 0, we find that x = 
Therefore, we obtain 



VVv^ 2 



i. 



V 



2 4 (a + l)V. 2 + l 



(iP 2 - l) 4 



(137) 



and i? 3 ( Q+1 )/ 4 = y/vp 2 -1/ip with ip > 1. Since w +oo 
when R — > 1, the tachyon field is not defined for all times. 

(v) For n = —2 and k < 0, we find that x = 
ip I \J\ — ip 2 . Therefore, we obtain 



*VpZ 



3c 2 
8nG 



1/2 



VT 



:E 3(l+a)/2 5 (M3) 



where i? = a/a*. Since p = p*/R^ 1+a ^ = 1/[6 tt(?(1 + 
a) 2 t 2 ], the scalar field evolves with time as <p = v 7 ! + a 
We note that F(</>) = for a = 0. On the other hand, 
there is no tachyon field corresponding to the vacuum 
equation of state p = —pc 2 (i.e. a = —1) used to describe 
inflation. By contrast, the equation of state (TfTj|) admits 
a tachyon field representation in the models with k < 
and a < exhibiting a phase of inflation. This may be an 
advantage of the polytropic equation of state (ITUl) with 
respect to the linear equation of state p = —pc 2 . 



V 2 



2 4 
~P*c 



(i - ip 2 Y 



(138) 



and i?3(a+i)/4 = y / l r^ 2 /^, with < V < 1- K a < 
0, w is between —1 and 0, so the tachyon field is well- 
defined. In particular, for a = 0, we have 



V(iP) 



p*c- 



(1- .02)3/2' 



(139) 



IX. CONCLUSION 

In this paper, and in the previous one, we have carried 
out an exhaustive study of the generalized equation of 
state (fTU)) . considering all the possible cases. We have 
obtained the following results (we assume here a ^ — 1): 

(i) For n > and k < 0, the universe undergoes an 
early inflation. It starts from t = — oo with a vanishing 
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radius and a finite density p ma x- Its radius increases with 
time while its density decreases. The universe exists at 
any time in the past and is non-singular. 

(ii) For n > and k > 0, the universe undergoes a new 
type of primordial singularity. It starts at t = with a 
finite radius and an infinite density. For t > 0, its radius 
increases while its density decreases. 

(iii) For n < and k < 0, the universe undergoes a late 
inflation. Its radius increases to +oo as t — > +oo while its 
density decreases to a constant value p m in ■ The universe 
exists at any time in the future and is non-singular. 

(iv) For n < and k > 0, the universe undergoes a 
future peculiarity. Its radius increases to a maximum 
value while its density decreases to zero (the universe 
"disappears"). For n < — 2, this peculiarity is reached 
in infinite time. For n > —2, this peculiarity is reached 
in a finite time (and there is a future singularity when 
n > — 1 due to the divergence of the pressure). Then, 
the radius decreases to zero while the density increases to 
+00. These phases of expansion and contraction continue 
periodically (cyclic universe). 

The generalized equation of state ([TO]) can be viewed as 
a "mixture" of a linear equation of state describing a clas- 
sical universe filled with radiation (a = 1/3) or matter 
(a = 0), and a polytropic equation of state whose ori- 
gin remains to be understood (in Paper I, we have men- 
tioned the connection with Bose-Einstein condensates, 
but other possibilities may be contemplated). Positive 
indices describe the early universe (the polytropic compo- 
nent dominates the linear component when the density is 
high) and negative indices describe the late universe (the 
polytropic component dominates the linear component 
when the density is low). A positive polytropic pressure 
(k > 0) leads to singular or peculiar models, while a neg- 
ative polytropic pressure (k < 0) leads to non-singular 
models that exist for all times. 

For k < 0, the generalized equation of state (|TU)) im- 
plies the existence of upper and lower bounds on the den- 
sity. It is natural to identify the upper bound with the 
Planck density p max — pp (quantum mechanics) and the 
lower bound with the cosmological density p m i n = p\ 
(general relativity). By taking n = 1 in the early uni- 
verse and n = — 1 in the late universe, we have obtained 
a non-singular model that is consistent with the known 
properties of our universe. This model improves the stan- 
dard ACDM model by removing the primordial singular- 
ity. An attractive feature of this model is its simplicity 
and the "symmetry" that it reveals between the past and 
the future. Furthermore, this model admits a scalar field 
representation based on a quintessence field (in the early 
and late universe) or a tachyon field (in the late universe). 

By using only the principle of "simplicity" , and reject- 
ing singular or peculiar solutions, we have obtained in 
a natural manner a simple model of universe. It cor- 
responds to an "aioniotic" universe that exists eternally 
in past and future and whose early and late evolutions 
present some striking symmetry. The remarkable point 
is that this model, which is obtained in a purely the- 



oretical manner independent of observations, turns out 
to be fully consistent with the known structure of the 
real universe. A theoretical challenge would now be to 
justify a quadratic equation of state p = —Ap 2 c 2 /(3pp) 
(n = 1) in the early universe and a constant pressure 
p = —pAC 2 (n = —1) in the late universe. It may also 
be interesting to study an equation of state of the form 
p/c 2 = — (a + l)p /pp + ap — (a + l)p\ exhibiting two 
phases of inflation [22j. For a = 1/3, it provides a uni- 
fication of pre-radiation, radiation, and dark energy (see 
Appendix |C|) . 

Our study suggests that the Planck density pp and 
the cosmological density p\ represent fundamental up- 
per and lower bounds. These bounds are responsible for 
a phase of early and late inflation. It is oftentimes ar- 
gued that the dark matter density p\ should be identi- 
fied with the vacuum energy. Since the vacuum energy 
is of the order of the Planck density pp, which is 10 122 
times larger than the cosmological density, this leads to 
the so-called cosmological constant problem [2(|. Actu- 
ally, if the Planck density and the cosmological density 
represent fundamental upper and lower bounds, it is not 
surprising that they differ by about 122 orders of magni- 
tudes. Therefore, the origin of the cosmological density 
should not be sought in quantum mechanics but in a new 
theory of cosmophysics based on general relativity. 

We have also studied singular models of universe cor- 
responding to a positive polytropic pressure (k < 0). Al- 
though these models do not appear to be selected by 
nature, they are interesting on a mathematical point of 
view. For n > 0, we have obtained a model exhibiting 
a primordial singularity at t = at which the radius 
has a finite value while the density and the tempera- 
ture are infinite. For n < 0, we have obtained mod- 
els exhibiting future peculiarities in which the density 
of the universe vanishes while its radius reaches a maxi- 
mum value. This peculiarity may occur in infinite time 
(n < —2), or periodically in time (n > —2). In partic- 
ular, we have obtained a simple analytical solution of a 
cyclic universe (anti-ACDM model) in which the density 
disappears and re-appears periodically. Since this solu- 
tion is always decelerating (in the phase of expansion), 
this peculiar model is not compatible with the known 
properties of our universe. 

There remains a type of solutions that we have not de- 
scribed so far: The case where the density increases as the 
universe expands (w < —1). These a priori un-natural 
solutions which violate the null dominant energ y c ondi- 
tion are associated with phantom scalar fields }23l |24| . 
They are studied in Paper III in which we construct 
models of "phantom universes" . Actually, there is a rich 
recent literature [f| [2f| on these solutions since obser- 
vations do not exclude the possibility that we live in a 
phantom universe [26|]. These models usually predict a 
future singularity at which the radius and the density of 
the universe become infinite in a finite time. This would 
lead to the death of the universe in a Big Rip or a Cosmic 
Doomsday [2~4j . 
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The main feature of the generalized equation of state 
(|10p is its poly tropic component p = kp 1 c 2 . A poly tropic 
equation of state occurs in many situations of astrophys- 
ical interest [![. For example, a poly tropic equation of 
state with index 7 = 5/3 describes adiabatic processes in 
main sequence stars [271 ] . Compact stars such as white 
dwarfs, neutron stars and BEC stars are also described by 
a polytropic equation of state. Classical white dwarf stars 
correspond to a polytropic index n = 3/2 and relativistic 
white dwarf stars correspond to a polytropic index n = 3 
[I?], [28| . Newtonian and semi- relativistic BEC stars with 
a self-interaction correspond to a polytropic index n = 1 
[29l | . From a mathematical point of view, polytropic stars 
represent a particular class of steady states of the Euler- 
Poisson system. Similarly, stellar polytropes represent a 
particular class of steady states of the Vlasov-Poisson sys- 
tem [l| . Polytropic distributions also appear in statistical 
physics, in relation with generalized thermodynamics (30j 
and nonlinear Fokker-Planck equations [3l| . In previous 
works, we have studied polytropic equations of state in 
several situations of astrophysical [321 [33| , physical [34| , 
and biological (3f| interest. It was therefore natural to 
consider the application of polytropic equations of state 
in cosmology also. 



reaches its maximum at 
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The evolution of the density and temperature as a func- 
tion of the scale factor is plotted in Fig. [17] 
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Appendix A: Equation of state p = (— p + kp~')c 2 with 
n < and k > 



FIG. 17. Evolution of the density and temperature as a func- 
tion of the scale factor. We have taken n = —1/2. 



In this Appendix, we specifically study the equation of 
state (|10l) with a = — 1, n < 0, and k > 0, namely 



p = {—p + kp J )c 2 . 



(Al) 



It generalizes the equation of state p = —pc 2 of the vac- 
uum energy. For the equation of state (|A1[) . the continu- 
ity equation ([7]) can be integrated into 



P = 



\n(a*/a) n ' 



(A2) 



where = (|n|/3fc) ra and a* is a constant of integration. 
The density is defined for a < a*. When a — > 0, p — > +00 
and p — > —00. When a — > a*, p — > 0. In the same limit, 
p — > for n < — 1, p tends to a finite value for 11 = —1, 
and p — > +00 for n > — 1. 

The thcrmodynamical equation ([5]) can be integrated 
into 



T = T» 



(n+l)/n 



(A3) 



where is a constant of integration. Combined with 
Eq. ([A"2|) . we obtain 



T 



ln(a*/a) n+1 Va 



(A4) 



When a -S- 0, T ->• 0. When a — s- a*, T — s- for n < -1 
and T — > +00 for n > — 1. For n < — 1, the temperature 



The equation of state can be written as p = wpc 2 with 

l/n 



M ( p 

3 VP* 

The pressure vanishes (11; = 0) at 



/0* 



a* 



(A7) 



(A8) 



When a — > 0, w — > —1; when a — > a*, w — > +00. The 
pressure is negative when a < a w and positive when a w < 
a < a*. 

The deceleration parameter is given by Eqs. (1-77) and 
(1-78). Together with Eq. (|AT| . we obtain 



(A9) 



The curve a(t) presents an inflexion point (a = q = 0) at 

(A10) 



p* 



— = e™ /2 . 
a* 



When a — ?► 0, q —> — 1; when a — >• a*, g — >• +co. The 
universe is accelerating when a < a c and decelerating 
when a c < a < a*. 

Finally, the velocity of sound is given by 



■i = - 1 - 



n+l f p 



l/n 



(All) 
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We have to distinguish two cases. For n > —1, the 
velocity of sound is always imaginary. When a — > 0, 



-1: when a 



a*, (c s /c) 2 



For 



n < — 1, the velocity of sound vanishes at the point (|A5[1 - 
(|A6|) where the temperature is maximum. At that point, 
the pressure is maximum with value 



Pe 

P*c 2 



3" 



(A12) 



When a — > 0, (c s /c) 2 — > — 1; when a — > a*, (c s /c) 2 — > 
+oo. The velocity of sound is imaginary when a < a e 
and real when a e < a < a*. On the other hand, the 
velocity of sound is equal to the speed of light at 



Ps 

p* 



G 



1 



a* 



(n+l)/6 



(A13) 



The velocity of sound is smaller than the speed of light 
when a < a s and larger when a s < a < a* . 

The evolution of w, q, and (c s /c) 2 as a function of the 
scale factor a is represented in Fig. [T51 




FIG. 18. Evolution of to, q, and (c s /c) 2 as a function of the 
scale factor a. We have taken n = —1/2. 



Setting i? = a/a*, the Friedmann equation © can be 
written 



7? 



e/ii? 



{-\nR) n / 2 ' 



(A14) 



where A" = (SnGp^/S) 1 / 2 and e = ±1 as explained pre- 
viously. We must distinguish three cases. 
For n < —2, 



R(t) = e-^™) 



-2/(|n|-2) 



-A't 
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(A16) 



The universe starts at t = with a vanishing radius, 
an infinite density, an infinite pressure, and a vanishing 
temperature. The universe expands and its density de- 
creases. It reaches its maximum radius R = 1 in infinite 
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FIG. 19. Evolution of the scale factor with time for n < —2, 
n = —2, and n > —2 (specifically n — —3, n — —2 and 
ji = -1/2). 



time, while the density and the pressure tend to zero 
algebraically rapidly. 
For n = —2, 



R(t) 



P(f) _ j$2 e -2Kt 



(A17) 
(A18) 



The universe starts from t = — oo with a vanishing radius, 
an infinite density, an infinite pressure, and a vanishing 
temperature. The universe expands and its density de- 
creases. It reaches its maximum radius R = 1 in infinite 
time, while its density and pressure tend to zero expo- 
nentially rapidly. 
For n > —2, 



R(t) 



p* 



-K(U - t) 



2M/(2-M) 
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(A20) 



The universe starts from t = — oo with a vanishing radius, 
an infinite density, an infinite pressure, and a vanishing 
temperature. The universe expands and its density de- 
creases. When t = £*, the universe reaches its maximum 
radius R = 1 while its density vanishes (future pecu- 
liarity). The universe "disappears". When t > 4*, the 
universe contracts (this corresponds to the solution of 
Eq. (|A14[) with e = — 1) and its density increases. When 
t — > +oo, the universe has a vanishing radius and an infi- 
nite density. At t = t*, the pressure vanishes for n < — 1, 
is constant for n = —1, and is infinite for n > — 1. In this 
last case, there is a future singularity. 

For n < —1, the velocity of sound becomes larger than 
the speed of light when t > t s corresponding to Eq. 
(|A13[) . This marks the end of the physical universe. The 
solutions (|A15|) - (|A20[) may not have sense for t > tt = t s . 
The evolution of the scale factor with time is plotted in 

Fig. m 
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Appendix B: The Planck scales and the cosmological 

scales 



It is expected that quantum mechanics plays a funda- 
mental role in the early universe where the scale factor 
is small. From the Planck constant K, the speed of light 
c, and the constant of gravity G, we can construct char- 
acteristic scales, known as the Planck scales. These are 
the Planck density 



Pp 

the Planck time 



G 2 h 



hG 

5 



5.1610 9y g/m 



the Planck length 
lp = 

the Planck mass 
M P 



Gh 



1/2 



1/2 



1/2 



5.3910" 44 s, 



1.6210 _35 m 



2.1710" 5 g, 



and the Planck temperature 



Mpc 2 

T P = — — = 1.42 10 32 K. 



(Bl) 



(B2) 



(B3) 



(B4) 



(B5) 



Mp/l 



3 

Pi 



We have the obvious relations lp — dp, pp 
tp = (Gpp) -1 / 2 . We have given in Paper I an inter- 
pretation of these different scales. In particular, the 
Planck density is expected to represent a fundamental 
upper bound for the density leading to a phase of early 
inflation. 

Similarly, we expect that the cosmological constant A 
plays a fundamental role in the late universe where the 
scale factor is large. We shall regard the cosmological 
constant as a fundamental constant of nature that de- 
scribes the cosmophysics (hence the late universe) ex- 
actly like the Planck constant describes the microphysics 
(hence the early universe). If this interpretation is cor- 
rect, we should not seek the meaning of the dark energy 
in the vacuum energy, or in quantum mechanics. It would 
have another origin purely connected to general relativ- 
ity. In this point of view, the cosmological constant prob- 
lem is a false problem (see Sec. I VII E[) . However, the true 
nature of the cosmological constant remains to be found. 
It would be necessary to develop the counterpart of quan- 
tum mechanics to cosmophysics. For the moment, there 
is no such theory. However, using dimensional analysis, 
we can easily introduce cosmological scales that are the 
counterpart of the Planck scales. 

From observations, we know the value of the dark en- 
ergy density p\. It is usually written as pa = ^A.oPo 
where 0,a,o — 0.763 is the present fraction of dark en- 
ergy and po = 9.20 10 _24 g/m 3 is the present density of 



the universe determined by the Hubble constant Hq = 
2.2710~ 18 s _1 through the Friedmann equation ©. This 
yields pa = 7.02 10 _24 g/m 3 . From now on, pa will be 
called the cosmological density (the word "cosmological" 
in the late universe replacing the name "Planck" in the 
early universe). From the cosmological density pa, the 
speed of light c, and the constant of gravity G, we can 
construct a cosmological time t\ = ^/(Gpa) 1 ^ 2 , a cosmo- 
logical length I a = ct\, and a cosmological mass Ma = 
PaI\- We have p A = 0, A ,oPo, *A = (Stt/SCIa^^Ho 1 , 
and J A = (87r/3f2 A ,o) 1/2 ao- 

The dark energy density is related to the cosmologi- 
cal constant by pa = A/(8nG). From now on, we will 



consider that A = 1.1810 



-35, 



is a fundamental con- 



stant of physics. Therefore, the corresponding scales pa, 
tA, Ia--- are fundamental scales determined from A in the 
same sense that the Planck scales are determined from K. 
These cosmological scales can be written as the Planck 
scales by introducing the notation 



8ttc 5 
~AG 



(B6) 



Numerically, Ti A = 7.35 10 122 ft = 7.75 1 88 Js. This con- 
stant might play a fundamental role in the theory of cos- 
mophysics, just like the Planck constant h plays a fun- 
damental role in the theory of microphysics. This is, 
however, a pure speculation that is simply motivated by 
the apparent "symmetry" between the early universe and 
the late universe (in any case, Ha can be introduced as a 
convenient notation) . The theory of cosmophysics would 
then lead to characteristic scales that are the cosmologi- 
cal density 



A 



PA 



G 2 h A 8ttG 
the cosmological time 

the cosmological length 

1/2 



7.0210" 24 g/m 3 , (B7) 
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-m - 



A 



1/2 



1/2 



the cosmological mass 

1/2 



M A 



Hac 
G 



(%-kAG 2 ) 1 ' 2 
and the cosmological temperature 
M A c 2 , 



1.4610 i8 s, (B8) 



4.3810 26 m, (B9) 



5.90 1 56 g, (BIO) 



10 93 K. 



(Bll) 



We have given in this paper an interpretation of these dif- 
ferent scales (except for the temperature whose nature is 
unknown to us). In particular, the cosmological density 
is expected to represent a fundamental lower bound for 
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the density leading to a phase of late inflation. Actu- 
ally, these scales were first introduced by Einstein [36| in 
his static model with a cosmological constant. Although 
they now have a different meaning (since the universe 
is not static), the scaling is of course the same. These 
scales, coming from the fundamental constant Ka (a re- 
definition of A) should ultimately emerge from a theory 
of cosmophysics that remains to be constructed. The ap- 
parent symmetry between the early and the late universe, 
may be a source of inspiration in the construction of such 
a theory. 



Appendix C: Unification of pre-radiation, radiation, 
and dark energy 

We have seen in Paper I that the transition between the 
pre-radiation era and the radiation era can be described 
by an equation of state p/c 2 = —4p 2 /(3pp) + p/3. On 
the other hand, the transition between a linear equation 
of state era and the dark energy era can be described by 
an equation of state p/c 2 = ap — (a + l)p\ (see Sees. IIIII 
and IIV[) . In Sec. |V] we have taken a = to describe the 
transition between the matter era and the dark energy 
era. Alternatively, if we take a = 1/3, we describe the 
transition between the radiation era and the dark energy 
era. Finally, the quadratic equation of state 



V 



v c2 

3pp 



;PA C 



(CI) 



describes, in a unified manner, the pre-radiation era, the 
radiation era, and the dark energy era. A nice feature 
of this equation of state is that both the Planck den- 
sity (early universe) and the cosmological density (late 
universe) explicitly appear. Therefore, this equation of 
state reproduces both the early inflation and the late 
inflation. This may suggest that the pre-radiation, radi- 
ation, and dark energy are the manifestation of a unique 
form of "generalized radiation". When p — > pp, we get 
p ~ —ppc 2 (vacuum energy); when p — > pa, we get 
P ~ —p\c 2 (dark energy); when p\ <C p <C pp., we get 
p ~ pc 2 /3 (radiation). 

For the equation of state (|C1[) . the continuity equation 
J7]) can be integrated into 



Pp 



i + (a/ ai y 



PA, 



(C2) 



where a\ = 2.61 10 6 m (see Paper I). To obtain this ex- 
pression, we have used the fact that pp p\ so that 
p/c 2 + p ~ (4/3p P )(p - Pa)(pp - p). When a -> 0, 
p — > pp (early inflation); when a — > +oo, p — > p\ (late 
inflation); when p\ -c p <C pp, p ~ pp(a/ai)~ 4 (power- 
law evolution). We note that Eq. (|C2I) can be rewritten 
as 



QradPO 



(a/a ) 4 + (ai/a ) 4 



(C3) 



It leads to Eq. (JST) with Cl m = Q B + Q DM = 0. It 
is therefore necessary to add "by hand" the density of 
matter (|94|) and (|95| in the Fricdmann equation ([9]) . This 
suggests that matter on the one hand, and pre-radiation 
+ radiation + dark energy on the other hand should be 
treated as two different "species". This is at variance 
with usual models, including those considered in the main 
part of this paper, that try to unify matter and dark 
energy. 

For the equation of state (|Cip . the thermodynamical 
equation ([5]) can be integrated into 



T = T P 



1 r \ 1/4 / 

15\ I p PA 
\Pp Pp 



1/4 



PPJ 



<7 4 



(CM) 



where we have determined the constant of integration in 
order to recover the Stcfan-Bolzmann law in the radiation 
era pa <C p <C pp (see Paper I). Combined with Eq. 
(EU), we get 



T = T f 



1/4 



(a/a,) 7 



[Kai) 4 + ir 



(C5) 



This is exactly the same expression as in the pre-radiation 
+ radiation era (see Paper I) but it is also valid in the 
dark energy era. This may be a confirmation that dark 
energy is of the same nature as pre-radiation and radi- 
ation. When o < oi, T ~ T P (15/7r 2 ) 1 / 4 (a/ai) 7 (pre- 
radiation); when a~^> a\, T ~ Tp(15/n 2 ) 1 ^ 4 (ai/a) (radi- 
ation and dark energy). 

Simple analytical expressions of the Friedmann equa- 
tion ([5]) with the density-radius relation (|C2[) can be 
given in particular limits. In the pre-radiation era, the 
scale factor increases exponentially as in Eq. (1-119). In 
the pre-radiation + radiation era, the evolution of the 
scale factor is given by Eq. (1-128); in the pure radiation 
era, the scale factor increases algebraically as in Eq. (I- 
124); in the radiation + dark energy era, the evolution 
of the scale factor is given by Eq. ([2~4")l with a = 1/4, 
a* = (pp/p A ) 1/4 ai, and K = (Sir/S) 1 / 2 ^ 1 ; in the dark 
energy era, the scale factor increases exponentially as in 
Eq. (|48p . Actually, it is possible to solve the Friedmann 
equation © with the density-radius relation (|C2[) ex- 
actly. Introducing R = aja\ and A = pa/ pp <C 1, we 
obtain 



1/2 



VTTWi dR=(^\ t/tp + C, (C6) 



which can be integrated into 
1 



In 



1 + 2\R A + 2^ X(l + R 4 + XR S ) 



In 



2 + R 4 + 2-v/l + R 4 + XR & 



R 4 



tp 



C, 



where the constant C is determined such that a 
t = (see Paper I). 



(C7) 
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